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ABSTRACT 

We present the fi-band Tully-Fisher relation (TFR) of 60 late-type galaxies with redshifts 0. 1 - L The galaxies 
were selected from the FORS Deep Field with a limiting magnitude of /? = 23. Spatially resolved rotation curves 
were derived from spectra obtained with FORS2 at the VLT. High-mass galaxies with Vmax ;J 150km/s show little 
evolution, whereas the least massive systems in our sample are brighter by ~ 1 - 2 mag compared to their local 
counterparts. For the entire distant sample, the TFR slope is flatter than for local field galaxies (-5.77 ±0.45 
versus -7.92 ±0.18). Thus, we find evidence for evolution of the slope of the TFR with redshift on the 3 a level. 
This is still true when we subdivide the sample into three redshift bins. We speculate that the flatter tilt of our 
sample is caused by the evolution of luminosities and an additional population of blue galaxies at z > 0.2. The 
mass dependence of the TFR evolution also leads to variations for different galaxy types in magnitude-limited 
samples, suggesting that selection effects can account for the discrepant results of previous TFR studies on the 
luminosity evolution of late-type galaxies. 

Subject headings: galaxies; evolution — galaxies: kinematics and dynamics — galaxies: spiral 
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1. INTRODUCTION 

Tully and Fisher (1977) found a remarkable correlation be- 
tween the stellar content and the kinematics of disk galaxies 
linking tightly the luminosity to the maximum rotational veloc- 
ity of late-type galaxies. This "Tully-Fisher Relation" (TFR) 
can therefore be exploited to study both the dynamical evolu- 
tion and the star formation history of disk galaxies. Observed 
local TFR samples (e.g. Mathewson and Ford 1996) with their 
very small scatter already served as important constraints for 
galaxy evolution models (e.g. Eisenstein and Loeb 1996). But 
the redshift evolution of the TFR is an even more powerful tool 
to investigate galaxy evolution and cosmological structure for- 
mation. 

On the observational side, recent studies of spiral galaxy 
kinematics at intermediate redshift produced rather discrepant 
results mainly due to different sample selection methods. While 
Vogt et al. (1996) find only a modest increase in luminosity 
of AMb « 0.5 at (z) « 0.5, the results of Rix et al. (1997) 
and Simard and Pritchet (1998) imply a strong brightening of 
roughly 2 mag in absolute B-band luminosity. 

Theoretical simulations have to consider the temporal growth 
of galaxies in size and mass as well as the evolution of the 
star formation rates. Mao et al. (1998), for example, combine 
the dynamical evolution as traced by Cold Dark Matter (CDM) 
hierarchical-merging models with a simple scaling relation be- 
tween mass and luminosity. Steinmetz and Navarro (1999), in 
addition, introduce star formation phenomena via semianalyti- 
cal recipes. These kind of studies predict zeropoint offsets in lu- 



minosity for distant TFR samples depending both on observed 
wavelength and cosmology. 



2. SAMPLE AND OBSERVATIONS 

Our galaxy sample is based on the FORS Deep Field (FDF), 
a deep multi-color study iUBgRIzJK) of an w 6' x 6' field with 
limiting blue and visual magnitudes comparable to the HST 
Deep Fields (Appenzefler et al. 2000; Heidt et al. 2001). Photo- 
metric redshifts and spectrophotometric galaxy types were de- 
termined for an /-band selected catalogue of « 3800 objects 
(Bender et al. 2001). From this list, late-type galaxies were 
drawn primarily on the basis of their brightness with a limit- 
ing magnitude of R < 23"' (as determined by the SExtractor 
package, Bertin and Arnouts 1996). Additional criteria have 
been that the inclination was / > 40° and Zphot < 1 .2 so that any 
[Oil] 3727 emission-line doublet would fall onto our observed 
wavelength range. The final sample comprised 79 spirals of all 
Hubble types (Sa - Irr) with 0. 1 < z < 1 .0 peaking at z = 0.35 
and 0.75 (see Fig. 1). Their absolute B magnitudes cover the 
range -22.4 < Mb < -15.8, the scale lengths range from 0.6 kpc 
to 4.8 kpc (Hq = 75 km s"'Mpc"' and qo = assumed throughout 



the paper). 

' Based on observations collected at the European Southern Observatory, Cerro Paranal, Chile (ESO Nos. 65.0-0049 & 66.A-0547). 
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Fig. 1 . — Redshift histogram of the FDF sample. Overplotted is the 
Mb distribution. 

Spectroscopy was carried out with the Focal Reducer and 
Spectrograph 2 (F0RS2) mounted onto the Very Large Tele- 
scope (VLT) Unit 2. Using the grism 600R with the order 
separation filter GG435, the typical wavelength range for the 
case of a chip-centered slit was AA = 5250 -7450 A. With a 
slit width of 1" the dispersion was 1.08 A/pixel and the resolu- 
tion R « 1230. Observations were done in MOS (Multi Object 
Spectroscopy) mode, in which FORS2 provides 19 sUtlets of 
« 22" length in a vertical array. Targets for the different setups 
were selected in a way to minimize the misalignment between 
slit direction and position angle of the galaxy. The maximal de- 
viation was 15°. The exposures for each setup were split into 
3 X 3000 s, the mean DIMM seeing was 0.66" (FWHM) for a 
spatial scale of 0.2"/pixel. A detailed description of the reduc- 
tion process will be given in Bohm et al. (2002). 

3. DETERMINATION OF Vmax AND Mb 

To derive the rotation velocities of the galaxies as a func- 
tion of radius, Gaussian profiles were fitted to the visible emis- 
sion line(s). For the objects with high redshift (0.4 < z < 1.0), 
the [Oii]3727 doublet was used, while for objects with 0.1 < 
z < 0.5, the [O iii]5007 and/or lines were measured. To in- 
crease the S/N, a boxcar filter of 0.6" was applied, i.e. three 
neighbouring rows of the spectra were averaged before each 
Gaussian fit. 

In contrast to studies of local spirals, these "straightforward" 
measurements of Vrot(r) cannot be used to directly derive the 
maximum rotation velocity v,n;,x (Simard and Pritchet 1999). 
Because the size of the visible disks of spirals at z fa 0.5 (with 
typical apparent scale lengths of r^ 0.5") is comparable to 
the slit width of 1", the spectroscopy is an "integration" over 
the galaxies' intrinsic velocity fields. 

To tackle this problem, a synthetic rotation curve (RC) was 
generated for each galaxy, assuming a linearly rising rotation 
velocity in the innermost galaxy regions and a constant rota- 
tion Vmax at large radii (e.g. Courteau 1997). The radius at 
which this intrinsic RC flattens was computed from the scale 
length Td measured in the /-band, taking into account that its 
value depends on the rest frame wavelength, because the scale 
length of the young star population is larger than that of the 
older population (Ryder and Dopita 1994). This intrinsic one- 
dimensional RC was then used to generate the 2D velocity field 
taking into account inclination and position angle. The struc- 
tural parameters were derived from two-dimensional luminos- 
ity profile fits based on a coadded I-band image with an FWHM 
of 0.5". These fits considered the possible dependence of the 
Point Spread Function on CCD position. If a bulge was de- 
tectable, the fits included two components. To test our pro- 
cedure, we performed a "VLT simulation" based on the driz- 
zled HDF images by rebinning them to 0.2"/pixel and convolv- 



ing them to 0.5" FWHM. The disc inclinations we measured 
in these frames were then compared with the values given by 
Marleau and Simard (1998) based on the original images. For 
objects with B/T<0.4, the differences were limited to a few de- 
grees. Estimated errors on the inclinations range from 7° for 
i « 40° to 3° for i f« 80°, in some cases they may be as high 
as up to 12°. In the next step, the inclined velocity field was 
folded with the luminosity profile and convolved to account for 
the seeing during spectroscopy. A simulated RC was extracted 
from the 2D velocity field by integrating over the slit aperture. 
At last, we variied Vmax. the only remaining free parameter in 
these simulations, until the best reproduction of the observed 
RC was achieved. In that manner we were able to quantify the 
maximum rotation velocity of 60 spirals in our sample, with es- 
timated errors between 10 and 40 km/s, depending on galaxy 
size and RC quality. The remaining 19 objects showed either a 
solid-body rotation, highly disturbed RCs or no rotation at all. 
Examples of RCs are given in Fig. 2. 
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Fig. 2. — Examples of observed rotation curves with their ID-fit. 
Respective values of redshift, /obs magnitude, disk scale length (in 
arcsec) and Vmax (as derived with the 2D-model of the velocity field) 
are given. Errorbars denote the errors of the Gaussian fits to the line 
positions. 

The 5-band luminosities of the FDF galaxies were deter- 
mined in the following way: we took SExtractor's MAG_ISO 
values as the total magnitudes in the BgRI FORS2 filters, which 
were then transformed onto the Johnson-Cousins photometric 
system. One origin of uncertainty in the derivation of the rest- 
frame fi-band luminosity is the ^-correction (see e.g. Ziegler 
et al. 1999). For early- type spirals at z w 1.0, transformation 
from comoving to rest-frame B-band magnitude amounts to 
> 2 mag. To keep the ^-corrections small, we therefore used 
that filter which best matches the rest-frame B. Depending on 
the galaxies redshift, this was either the B-, g-, R- or /-band, 
leading to a typical correction of 0.5 mag. For the calculation 
of the respective ^-corrections, we used template spectra cre- 
ated by the chemically consistent evolutionary synthesis mod- 
els of MoUer et al. (2001). Observed galaxy colors of different 
Hubble types (e.g. Frei and Gunn (1994)) are reproduced with 
these spectra to within 0.15 mag over the whole covered red- 
shift range. Corrections for inclination-dependent intrinsic ab- 
sorption due to dust were appUed following Tully and Fouque 
(1985). 
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4. THE TULLY-FISHER RELATION AT Z « 0.5 

For the comparison to local data, we chose the large sample 
of 1 196 field Sbc and Sc galaxies of Haynes et al. (1999). To 
be able to better compare our results with previous studies, we 
construct the B-band IFK in this letter (for the /-band and NIR 
TFR see Bohm et al. 2002). Therefore, the observed / magni- 
tudes given by Haynes et al. were transformed to B using the 
type-dependent (B-I) colors of Frei and Gunn (1994) and cor- 
rected for intrinsic absorption the same way as our sample. This 
transformation was tested on our sample (for which we both 
have B and / magnitudes), showing that on average the accu- 
racy is within 0.1 mag independent of galaxy type. The result- 
ing absolute magnitudes cover -22.5 < Mb < -15.1. For use in 
the TFR, we further corrected the derived Mb magnitudes for 
morphological and incompleteness bias following Giovanelli 
et al. (1997). The authors pointed out that an incompleteness 
bias arises from the differences between the luminosity function 
of a magnitude limited sample and the Schechter luminosity 
function. At low rotation velocities (low absolute magnitudes), 
galaxies falling below the detection limit do not contribute to 
the observed TFR, resulting in an underestimated slope. Sim- 
ulating this, a synthetic, unbiased TFR is transformed into the 
observed relation using the measured values of Vmax, the ob- 
served scatter (as a function of v^ax) and the completeness limit. 
The effect can then be counter-balanced by decreasing the ob- 
served luminosities accordingly. 

The magnitudes of the distant galaxies were de-biased with 
the same procedure adopting the parameterization of the scatter 
given by Giovanelli et al. but taking into account that the scatter 
in the FDF sample is on the mean 1 .5 times larger. To account 
for the variation in the luminosity function with redshift, we 
subdivided our full FDF sample into three z bins, for which we 
constructed the observed luminosity functions separately. At 
the faint end of our sample, the correction for incompleteness 
amounts to maximal +1.4 mag. The morphological bias, as a 
similar effect, originates in the different luminosity functions 
for different morphological classes and leads to a slightly in- 
creased luminosity for early types (~ -0.3... -0.1 mag). If not 
accounted for, the combination of these two biases leads to an 
underestimate of both slope and scatter of the TFR. We sub- 
divided the 60 galaxies into 3 classes corresponding to Sa/Sb, 
Sc and Sd/Irr. This was done on the basis of spectrophotomet- 
ric types from our photometric redshift catalog and by compar- 
ing our spectra to the catalog of template spectra of Kennicutt 
(1992). There are 10 Sa/b, 32 Sc and 18 Sd/Irr galaxies in our 
sample. 

We present the distant B-band TFR in Fig. 3 in comparison 
to the local TFR. The FDF spirals follow a flatter relation than 
their local counterparts. While the most massive galaxies are 
indistinguishable on average from the local sample, the less 
massive galaxies show large deviations from the local fit. For 
the full FDF sample we find a slope of -5.77 ± 0.45 by repeat- 
ing bisector fits (e.g. Ziegler et al. 2001) to the data centered 
on the median value of the logVmax-distribution ((logVmax) = 
2.08) 100 times using a bootstrap method. For the local sam- 
ple ((log Vmax) = 2.182) the same method yields a slope of 
-7.92 ±0.18. The last value is consistent with slopes derived 
for other local field samples. We derive, for example, a slope of 
-7.40 for the BBFN/RC3 sample (not de-biased; Burstein et al. 
1997), and -7.50 for the Pierce and Tully (1992) sample. Sakai 
et al. (2000) find a slope of -8.07 ± 0.72 in the B-band using 
HST-based calibrators. Thus, the slopes of the local samples 



and the FDF sample differ on the 3 a level. This is confirmed 
by a 2D Kolmogorov-Smirnov test, which limits the probabil- 
ity that our sample originates from the same distribution as the 
Haynes et al. sample to 0.001. We also performed a simple 
Monte-Carlo simulation extracting randomly 60 galaxies from 
the local sample. Neither of 200 iterations resulted in a TFR 
slope as flat as for the FDF galaxies. If we subdivide our sam- 
ple according to galaxy type (Sa/b, Sc and Sd/Irr) or redshift, 
the resulting slopes and zero-points are compatible with each 
other within their la errors (see Table 1). 




log v^^^ (km s-i) 

Fig. 3. — The S-band Tully-Fishcr relation of 60 field spiral and 
irregular galaxies at intermediate redshifts. Big symbols denote our 
spectrophotometric classification: squares Sa/Sb, circles Sc, triangles 
Sd/Irr. Open symbols correspond to "observed" data, filled symbols to 
de-biased data. Solid (dotted) lines represent the best (± 1 cr) fits of the 
100 bootstrap bisector fits to the de-biased sample. The dashed lines 
are the fits to the "observed" data. The steeper lines show similar fits 
to the large sample of 1200 local field spirals (dots) from Haynes et al. 
(1999). There is a clear change of slope (on the 3 a level) between the 
local and distant TFR. Two typical errorbars are given. 



Table 1 
tfr slopes with 1 a errors 



sample 


N 


Mb vs. logv 


logv vs. Mb 


distant sample: 








FDF (not de-biased) 


60 


-4.99 ±0.48 


-0.202 ±0.016 


FDF 


60 


-5.77 ±0.45 


-0.175 ±0.014 


FDF z < 0.3 


11 


-5.46 ±0.68 


-0.195 ±0.054 


FDF 0.3 < z < 0.5 


29 


-4.62 ±0.48 


-0.220 ±0.026 


FDF z > 0.5 


20 


-5.59 ±0.61 


-0.174±0.014 


local samples: 








Haynes et al. (not de-biased) 


1296 


-7.52 ±0.18 


-0.133 ±0.003 


Haynes etal. (1999) 


1296 


-7.92 ±0.18 


-0.126 ±0.003 


BBFN (not de-biased) 


467 


-7.40 + 0.21 


-0.135 + 0.004 


Sakai el al. (2(K)()) 


21 


-S.()7 1 0,72 





Can the observed change of slope be caused by systematic 
errors? Any overestimate of / would both lead to an underes- 
timate of Vmax and too large a luminosity after correcting for 
intrinsic absorption. But to introduce a signiflcant effect on the 
slope, the errors on / would have to be large: assuming, e.g., an 
inclination of 60° instead of 40° for an object, the shift in the 
TFR with a slope of -6 would amount to -0.93 mag. Such mea- 
surement errors are unlikely, as our incUnations do not show 
any dependence on rotation velocity nor apparent scale length 
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nor type. For that reason, we are also confident that morpho- 
logically irregular objects in our sample are unlikely to have 
introduced a significant bias in the inclinations. 

5. DISCUSSION 

Our intermediate-z TFR indicates a change in the slope with 
look-back time. Most previous studies of the kinematics of 
spirals in this regime were restricted to small data sets (A^ < 20) 
and had to assume a constant slope. In the ongoing study of 
Vogt (2000) of ^ 100 distant disk galaxies, it was also assumed 
that the slope of the local TFR of Pierce and TuUy (1992) holds 
vaUd at all redshifts. 

CDM-based simulations predict a slight brightening of the 
B-band TFR up to z = 1 (e.g. Steinmetz and Navarro 1999). In 
this model, the decrease in mass is more than balanced by an 
increase in luminosity due to younger stellar populations. Any 
mass-dependent evolution of the star formation rate or other TF 
parameters that would cause an evolving TFR slope are yet not 
implemented in these models. 

A change in the TFR slope as we observed makes it impos- 
sible to give a single number for the luminosity evolution of 
distant galaxies since it depends on the galaxy mass. The FDF 
galaxies with v^.^^ > 150km/s fall into the same magnitude re- 
gion as those of the local sample and show no significant lumi- 
nosity evolution, whereas the less massive distant galaxies are 
brighter by ~ 1-2 mag. 

A mass dependence of the TFR evolution also introduces 
a dependence on galaxy type since the various sub-classes of 
disk galaxies have different distributions in mass (logVmax) for 
a magnitude-Umited sample. Thus, a natural explanation of the 
discrepant results in luminosity evolution of previous observa- 
tional studies is suggested. For example, studies based on sam- 
ples with galaxies selected to be blue (hke Rix et al. 1997) or 
showing strong emission lines (like Simard and Pritchet 1998) 
are biased towards later types of spirals (de Vaucouleurs 1961; 
Kennicutt and Kent 1983) resulting in a large luminosity in- 
crease. On the other hand, studies using mainly galaxies with 
large disks (hke Vogt et al. 1996) may contain a larger fraction 



of early and luminous spirals leading to a more modest lumi- 
nosity evolution. 

At last, we speculate on the possible origin of the flatter tilt 
in all three redshift bins (z < 0.3,0.3 < z < 0.5,z > 0.5). A 
plausible scenario would be that we not only trace a pure lu- 
minosity evolution but also observe a population of blue dwarf 
galaxies with Mb > -20 underrepresented in the local Universe. 
This would be consistent with blue number counts that suggest 
a higher number density of bursting dwarfs at z < 0.5 (e.g. El- 
lis 1997). These objects may have too low a rotation for their 
luminosity as observed by Rix et al. (1997). On the other hand, 
the star formation rates of all FDF galaxies (0.2- ISM© /yr) 
fall well below the local regime of starburst galaxies (Kennicutt 
1998). If we exclude the low-mass objects with v^ax < 80km/s 
and Mb > -18 (A' = 16) from the local and distant sample, we 
derive slopes of -7.29 ±0.15 and -5.08 ±0.50, respectively, 
confirming the change of tilt with 3 a confidence. 

6. SUMMARY 

We have presented the B-band TuUy-Fisher relation of 60 
disk galaxies at intermediate redshifts. The galaxies comprise 
a magnitude-limited sample within the FORS Deep Field. We 
observe a significant (3(t) change of slope in comparison to 
the local sample. The larger FDF spirals (vmax ^ 150km/s) fall 
into a region of the logVmax - luminosity diagram consistent 
with modest or no luminosity evolution. On the contrary, the 
smaller distant galaxies are offset from the local TFR by order 
~ 1 — 2mag. 
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Foundation (1/76520), the Deutsche Forschungsgemeinschaft 
(SFB375, SFB439) and the German Federal Ministry for Ed- 
ucation and Science (ID-Nos. 05 AV9MGA7, 05 AV9WM1/2, 
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